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Abstract 

The possibility of coherent radiation by molecular magnets is investigated. It is 
shown that to realize the coherent radiation, it is necessary to couple the considered 
sample to a resonant electric circuit. A theory for describing this phenomenon is de- 
veloped, based on a realistic microscopic Hamiltonian, including the Zeeman terms, 
single-site anisotropy, and dipole interactions. The role of hyperfine interactions be- 
tween molecular and nuclear spins is studied. Numerical solutions of the spin evolution 
equations are presented. 



PACS: 

07.57.Hm - Infrared, microwave, and radiowave sources 
75.50.Xx - Molecular magnets 

76.20. +q - General theory of magnetic resonances and relaxation 



Molecular magnets is a relatively new class of materials possessing several interesting 
properties (see e.g. reviews [1-3]). One of the important features is that the molecules, 
forming the magnet, can have rather different values of spin, between 1/2 and more than 10. 
At low temperatures, the molecular magnets are bistable and exhibit hysteresis [4] and tun- 
neling between up and down orientations [5,6]. Their potential use for quantum computation 
has been proposed [7], and experiments with millimeter- wave radiation have been carried out 
showing that the relaxation rate for magnetization reversal and the energy-level populations 
can be controlled [8-12]. The magnetization reversal, caused by a linearly varying in time 
external magnetic field, with a fast magnetization avalanche [4], is due to the Landau-Zener 
tunneling [13-16]. The latter is a quantum effect, involving no coherence between different 
molecules. An attempt has been made [17] to check if the quantum tunneling can be ac- 
companied by coherent electromagnetic radiation. However, the most recent and accurate 
experiment [18] detected no significant radiation. 

In the present paper, we suggest a principally different setup for realizing coherent ra- 
diation from molecular magnets. The key difference is that a magnet has to be coupled to 
a resonant electric circuit. Such a setup has proved to be effective in achieving coherent 
radiation by nuclear spins, which was described theoretically [19-21] as well as performed 
experimentally (see reviews [22-24]). Molecular magnets, however, have several properties 
essentially distinguishing them from nuclear magnets. The most important distinctions are 
the possibility of possessing a high spin, the existence of level splitting, the presence of a 
single-site anisotropy that can be rather large, and the necessity of taking account of the 
line-width narrowing effect due to a strong spin polarization. The aim of this letter is to 
present a generalized theory of coherent spin radiation, including these principal features, 
typical of molecular magnets, and to give numerical solutions of the derived generalized 
equations. 

Our consideration is based on the macroscopic Hamiltonian 

# = E^ + ^E^> (!) 

i i^j 

describing N molecules, each having spin S. The first term in Eq. (1) contains the single-spin 
Hamiltonian 

Hi = -fioB ■ S t - D{S*) 2 , (2) 

in which fiQ = h'ys, with 75 being the gyromagnetic ratio, B is the total magnetic field acting 
on spin Sj, and D is the anisotropy parameter. The total magnetic field 

B = B e z + (h + hi cos uot + H)e x (3) 

consists of a longitudinal magnetic field B , transverse constant field h , modelling the level 
splitting, an alternating field hi cos out, and the resonator field H. The latter is defined by 
the Kirchhoff equation 

— + 2 7 # + u 2 j Q H(t') dt' = 7e/ - 47T77 , (4) 

in which 7 is the circuit damping, u is the circuit natural frequency, e/ = h 2 cos out is an 
electromotive force, rj is a filling factor, and m x = (fi /V) J2i < Sf > is the average transverse 
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magnetization. The coupling of a molecular magnet to a resonant electric circuit through 
the Kirchhoff equation (4) is the fundamental point in the proposed setup. The second term 
in Hamiltonian (1) corresponds to dipole spin interactions, with 

Hv^DgSfSf, (5) 

a/3 

,a(3 



where D"- is the dipolar tensor. 



The necessary conditions for the occurrence of coherent spin motion are the existence of a 
well-defined frequency of their rotation oo s and the formation of spin packets of a characteris- 
tic size L s , such that kL s <C 1, where k = uj s /c. Then the arising coherent radiation happens 
at the frequency uj s . If the radiation wavelength is larger than the system linear size L, then 
the spin-packet length L s , just coincides with L. This, however, is not compulsory, and L s 
can be much shorter than L. But the condition kL s <C 1 is necessary for spin packets to be 
formed. Writing down the Heisenberg equations of motion for the spin operators, we shall 
average them over the spin degrees of freedom, aiming at obtaining the evolution equations 
for the following averages: The transition function 



U 



E <sr>, (6) 



the coherence intensity 

w 

and the spin polarization 



S*N S (N S - 1) g < StS > > ' (7) 



N s 



' s mZ <S!> ' (8) 

where N s = pL? s is the number of spins in a spin packet, p = N/V is the average spin density, 
and Sf are the ladder spin operators. 

In each spin system, there are dipole spin interactions, represented by Eq. (5). These 
interactions must be taken into account when considering any collective properties of the 
spin system. The dipole interactions is the major cause for the appearance of the transverse 
relaxation rate. A peculiarity of molecular magnets is the possibility of having very high 
longitudinal spin polarizations. This requires to include in the effective transverse relaxation 
the effect of line- width narrowing [25], which results in the effective relaxation rate 

r 2 = (i- s 2 ) 72 + 7 2 *, (9) 



where 72 = n^pp^^J S(S + l)/h is the spin dephasing rate, with no being the number of 
nearest neighbours, and 72 is an inhomogeneous broadening [25,26]. 

From the evolution equations for the spin variables (6) to (8), it follows that the frequency 
of spin rotation is 

U s = UJ - UJ D S , (10) 

in which uo = —poB /h is the Zeeman frequency and u>d = (2S — 1)D /% is the anisotropy 
frequency. The second term in Eq. (10) depends on time through s = s(t). Therefore to keep 
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the spin frequency (10) in resonance with the resonator natural frequency u, one has to apply 
a sufficiently strong external magnetic field B , such that \u /u D \ ^> 1, and the resonance 
condition |A/u;| <C 1, (A = uj — \uj s \) could be realized. In some molecular magnets, the 
single-site anisotropy is rather strong, so that ujd can reach the values of 10 12 s _1 . Then, 
to make much larger then up, one has to invoke magnetic fields B ~ 10 — 100 T. 
Fortunately, many molecular magnets possess much weaker anisotropics, for which B ~ 1 
T would be quite sufficient. In addition, achieving nowadays very high magnetic fields is 
not an insuperable obstacle. Among available sources [27], there are those allowing to reach 
the stationary magnetic fields of 45 T (National High Magnetic Field Laboratory, USA) 
and the pulsed fields up to 600 T (University of Tokyo, Japan). The pulsed fields can be 
supported during the time 10~ 2 s up to several seconds, which is perfectly sufficient for 
realizing coherent radiation in a molecular magnet. 

In what follows, we assume that there is an external constant magnetic field B , for 
which the resonance condition |A/a>| <C 1 is valid, that is the resonant circuit can be tuned 
to the resonance with the spin frequency (10). As is mentioned above, to achieve coherent 
radiation, one has to have a well-defined frequency u s , which requires that the magnitude 
of the latter be much larger than any attenuation rate. If so, the system possesses two 
different time scales, which allows us to employ the scale-separation approach [19-23], whose 
mathematical foundation is based on the averaging techniques [28] . The transverse function 
(6) is fast, as compared to the slow variables (7) and (8). The latter are to be treated as 
quasi-invariants, when solving the equation for function (6). 

To represent the final formulas in a nice way, we need to introduce several notations. 
The spin-resonator coupling is characterized by the dimensionless coupling parameter 

9= , 2 , A2 v , (11) 
72 (7 2 + A 2 ) 

in which 7 = Tcripji^S '/% is the natural spin width. The action of the resonator feedback 
field comes through the temporal coupling function 

a = g l2 (l - e"*) . (12) 

The effective transverse attenuation is given by the collective width 

r = r 2 + 73 -as, (13) 

in which 73 is the inhomogeneous dynamic broadening due to local dipole fluctuations [19- 
23]. 

For the transverse function (6), we find 



uj s -iT A + iT 



Uo + uj s -tT A + iT 



where v = /i ho/h, v x = /i hi/2h, u 2 = //o^2/2^, and the function (3 = (i/ 2 /2)(l — e~ 7 *) 
describes the action of an electromotive force, if any. 
Let us define the effective attenuation 

r 3 = 73 + _ "o^ + W e _ rt (^ + ^) 2 r fl _ „x 
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Then the final evolution equations for the slow functions (7) and (8) can be represented in 
the form 



in which ( is a stationary spin polarization and Ti = 71 + 7* is the sum of the spin-lattice 
attenuation 71 and a pumping rate 7^, if the sample is subject to a nonresonant pumping 
procedure. 

The evolution equations (15) include all important attenuation rates that can influence 
spin dynamics. In general, there exists also one more relaxation parameter, the radiation rate 
7 r = 2fiQSk 3 N s /3h due to spin interactions through the common radiation field. However 
this rate, as has been noticed yet by Bloembergen [29], does not influence the dynamics of a 
macroscopic spin sample. This is evident from the ratio r ) r /l2 ~ 0.1 (kL s ) 3 <C 1, which shows 
that 7 r is negligible as compared to the line width formed by the dipole interactions. Such a 
situation is drastically different from atomic systems [22], in which both the dephasing rate 
72 as well as the collective radiation rate 7 r are caused by the same reason, by atomic inter- 
actions through the photon exchange, so that 7 r /72 ~ N c 1, with N c being the number 
of correlated atoms in a wave packet. But in spin systems, 7 r and 72 are due to different 
origins. As is thoroughly explained in review [24], in resonatorless spin systems interacting 
through dipole forces, the effect of pure spin superradiance is principally inacheivable. This 
is why, for organizing coherent radiation by a magnet, one has to couple the latter to a 
resonant circuit. 

Let us investigate the spin dynamics, described by Eqs. (15), for a molecular magnet, 
such as Mni 2 or Fe 8 , whose molecules possess spin S = 10. Then w D ~ 10 12 s _1 . At low 
temperatures, below the blocking temperature of about 1 K, the sample can be polarized, 
keeping its polarization for rather long times of order T\ ~ 10 5 — 10 7 s, hence 71 = Tf 1 is 
very small, 71 ~ 10~ 7 — 10~ 5 s -1 . The density of molecules in a magnet is p ~ 10 20 cm~ 3 . 
Dipole interactions are quite strong, yielding 72 ~ 10 10 s _1 . The tunnel splitting is small, 
with z/ ~ 10 4 s -1 , which is much smaller than uo- 

Assume that the sample is polarized, with the spin polarization up, Sq = s(0) > 0. Then 
it is placed in an external magnetic field B Q , such that the Zeeman frequency be positive, 
ujq > 0, hence /io-Bo < 0. Since fiQ = — 2/i#, where hb is the Bohr magneton, we have 
PbBo > 0. This means that the magnetic field is directed up, which implies that the magnet 
is prepared in a nonequilibrium state. The relaxation from this state is governed by Eqs. 
(15). We have analysed these equations numerically for a wide range of parameters. Our 
main concern has been to study collective effects without imposing transverse fields. So, we 
present here the results of our calculations, when v\ = (3 = 0, as a consequence of which 
r3 ~ 73- We assume that there is no pumping, so that 7* = 0, and r\ = 71 <C 72. The 
coupling parameter (11) is of the order of the resonator quality factor, g ~ Q. We have 
investigated the behaviour of the solutions to Eqs. (15) for the parameters 72, 73, 7, and 
g varied in a very wide diapason, up to three orders of magnitude, and for different initial 
conditions wq and sq- For a strong spin-resonator coupling g 1, all solutions, under fixed 
initial conditions, are qualitatively the same. Therefore, it is sufficient to demonstrate typical 
classes of solutions for a fixed set of the parameters 72, 73, 7, and g, but for distinct initial 
conditions. It is possible to distinguish three such qualitatively different classes of solutions, 
for each of which the relaxation time is much shorter than T 2 = 1/72- 




(15) 
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Figure 1 demonstrates the pure spin superradiance, starting from a high initial spin 
polarization s = 1, and without any initial coherence imposed on the sample, when w = 0. 
This means that superradiance develops as a purely self-organized process. Figure 2 shows 
the process of triggered spin superradiance, when the initial spin polarization s = 0.7 is yet 
high, an essential self-organization is present, but the relaxation is triggered by an initially 
induced coherence wo = 0.51. And Fig. 3 presents the effect of collective spin induction, in 
which the initial spin polarization s = 0.1 is low, playing not so important role, while the 
relaxation is induced by a high initial coherence wq = 0.99. In all the cases, the relaxation 
time, much shorter than T 2 , is due to collective effects caused by the resonator feedback field, 
with a strong spin-resonator coupling g ^> 1. If the latter is weak, g < 1, then collective 
coherent effects do not develop, and the relaxation time is of the order of T 2 . Recall that 
the initial coherence function w = \u \ 2 is nothing but the amplitude of the transverse spin 
squared. Therefore to create a nonzero Wq it is sufficient to get a nonzero transverse spin 
polarization uq, which can be obtained by imposing at the initial time a sufficiently large 
transverse pulsed field. 

In molecular magnets, in addition to dipolar fields, there also exist hyperfine fields due 
to interactions between molecular and nuclear spins. The influence of the latter can be 
taken into account by generalizing the Hamiltonian (1) to the form H = Hs + Hi + His, m 
which Hs is the same Hamiltonian (1) due to S-spins, the second term is the Hamiltonian 
of nuclear spins I, 

% i^j a/3 

where D°j} is a dipolar tensor for nuclear spin interactions, and the hyperfine interactions 
are described by the Hamiltonian 

i l^j a/3 

where A is the intensity of single-site hyperfine interactions and A?? is the dipolar tensor of 
hyperfine interactions at different sites. The following consideration can be done in the same 
way as earlier, by employing the scale separation approach [19-23]. Hyperfine interactions 
result in the appearance of an additional line broadening of the order of 10 7 — 10 8 s _1 , which 
is much smaller than the dipolar line width 72 ~ 10 10 s _1 . 

The most nontrivial effect caused by the presence of the hyperfine interactions is that the 
ensemble of nuclear spins, being also coupled to an electric circuit, serves as an additional 
resonator for molecular spins. This changes the effective coupling of molecular spins with 
the resonator circuit from Eq. (11) to the value 

= 77o^5 L pHiAsil \ 
72 (7 2 + A 2 ) \ p^ohuji ) 

where u>s = ujq — oops + (A/h)sjS and ujj = uj j + (A/h)sI, with uj j = —fijB /h being the 
nuclear Zeeman frequency, and sj is an average nuclear spin polarization defined analogously 
to Eq. (8). Estimates show that the change in the coupling parameter g, owing to the 
presence of nuclear spins, is very small, of the order of 10~ 5 . Thus, hyperfine interactions do 
not much influence the collective effects of molecular spins. 
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The intensity of radiation from each spin packet is I(t) = (2/3c 3 )filS 2 uj 4: N^w(t), which, 
depending on the parameters, can reach very high values of many Watts. This microwave 
radiation can be used for creating spin masers [3]. It is also feasible to organize the regime 
of punctuated spin superradiance [30], which can be employed for information processing. 
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Figure 1: Pure spin superradiance. Typical behaviour of the coherence intensity w(t) (solid 
line) and of the spin polarization s(t) (dashed line) as functions of time (in units of T 2 ) for 
the following attenuation parameters (in units of 72): 7^ = 1, 73 = 0.1, and 7 = 10. The 
spin-resonator coupling is g = 100. The initial conditions are w = and s = 1 
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Figure 2: Triggered spin superradiance. The coherence intensity w(t) (solid line) and spin 
polarization s(t) (dashed line) versus time (in units of T 2 ) for the same parameters as in Fig. 
1, but for the initial conditions w = 0.51 and s = 0.7 
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Figure 3: Collective spin induction. The coherence intensity w(t) (solid line) and spin 
polarization s(t) (dashed line) as functions of time (in units of T 2 ) for the same parameters 
as in Fig. 1, but for the initial conditions w = 0.99 and s — 0.1 
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